How tight is the Lieb-Oxford bound? 
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Density-functional theory requires ever better exchange-correlation (xc) functionals for the ever 
more precise description of many-body effects on electronic structure. Universal constraints on 
the xc energy are important ingredients in the construction of improved functionals. Here we 
investigate one such universal property of xc functionals: the Lieb-Oxford lower bound on the 



where C < Clo 



To this end, we 



exchange-correlation energy, S xc [n] > — C J d 
perform a survey of available exact or near-exact data on xc energies of atoms, ions, molecules, 
solids, and some model Hamiltonians (the electron liquid, Hooke's atom and the Hubbard model). 
All physically realistic density distributions investigated are consistent with the tighter limit C < 1. 
For large classes of systems one can obtain class-specific (but not fully universal) similar bounds. 
The Lieb-Oxford bound with Clo = 1-68 is a key ingredient in the construction of modern xc 
functionals, and a substantial change in the prefactor C will have consequences for the performance 
of these functionals. 



I. INTRODUCTION 

Any numerical calculation of the electronic struc- 
ture of matter that uses density-functional theory 
(DFT) employs an approximate exchange-correlation 
(xc) functional i 2 -! 3 - Further progress in DFT thus de- 
pends crucially on the development of ever better density 
functionals. A most important ingredient in this quest 
for better functionals is the small, but increasing, list of 
exact properties of and constraints on the universal xc 
functional 

In particular, this functional, E xc , is known to satisfy 
the following inequalities 4 ' 5 



> E xc [n] > B[n] > -C J dt 



r n 
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(1) 



where C 

,-i 



is a 



universal constant, B[n] = 
linXy—kO 7 -1 -E xc [n»y], and n 1 is the scaled density 
n 7 (r) = 7 3 n(7r). While the first inequality, providing 
an upper bound on E xc [n], is an immediate consequence 
of the variational principle, the second and third, 
providing lower bounds, are more complex. The second 
inequality, the Levy-Perdew boundf 5 - is based on scaling 
arguments. It contains E xc on the left and on the 
right, and thus provides a consistency test for any given 
approximation to E xc [n] . 

The third inequality is a remarkable result due to Lieb 
and Oxford,—, who established the form of the bound and 
obtained the value Clo = 1-68 as an upper limiti of the 
prefactor C. The present work is mostly concerned with 
this Lieb-Oxford bound, although a numerical compari- 
son with the Levy-Perdew bound will also be given. In 
terms of the local-density approximation (LDA) to the 
exchange energy, 
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the LO bound can also be written as 4 : 5 
E xc [n] > \LoE x DA [n\, 



(3) 



where Xlo = 1.354Cxo = 2.275. The analysis below is 
couched in terms of A. 

The Lieb-Oxford lower bound on the xc energy is one 
of not many exactly known properties of the univer- 
sal xc functional. Similarly to other such properties, is 
has been used as a constraint in the construction of ap- 
proximations to this functional! 4 " It is satisfied, e.g., by 
the LDA;£ the PBE generalized-gradient approximation 9 
(GGA) and the TPSS meta-GGA^ On the other hand, 
earlier GGAsii and semiempirical functionals containing 
fitting parameter s 12 ' 13 are not guaranteed to satisfy the 
bound for all possible densities. 

Note that Eq. ([3]) is a bound in the mathematical 
sense, i.e., -Ea; C [n] can never be more negative than 
\LoE x DA [n]. It is, however, not clear from the inequal- 
ity itself if Xlo is the smallest possible value of the pref- 
actor, i.e., if the bound can be tightened or not. Indeed, 
Chan and Handy 1 ^ have revisited the original calcula- 
tion of Lieb and Oxford, and obtained the slightly tighter 
bound X C H = 2.2149 (or C C h = 1.6358). 

Independently of the question whether the bound can 
be tightened mathematically, it is not clear if nature actu- 
ally makes use of the entire range of values of E xc allowed 
by the bound, and neither how distant specific classes of 
actual physical systems are from the mathematical max- 
imum. To put these issues in clearer focus, note that for 
any actual density n(r) one can, in principle, evaluate 
the density functionals E xc and E X D on this density, 
and calculate the ratio 



Ex DA [r 



= : X[r. 



(4) 



(2) 



which measures the weight of LDA exchange relative to 
the full exchange-correlation energy. The resulting value 



2 



of A[n] must be smaller than or equal to Xlo = 2.275, 
for any n, but it is not a priori clear by how much, and 
neither what the variations of A over different classes of 
systems are. 

Such information is not easy to obtain, since the def- 
inition (j4|) requires knowledge of the exact xc energy in 
the numerator and of the exact density n(r) in the nu- 
merator and the denominator. This knowledge is, in gen- 
eral, not available. There are, however, certain classes of 
systems for which near-exact xc energies and densities 
are available, e.g., from quantum Monte Carlo (QMC) or 
configuration interaction (CI) calculations. In this work 
we present a survey of avaliable such data for large and 
distinct classes of systems, and confront the results with 
the Lieb- Oxford bound. 

Section |TT] deals with real atoms, Sec. IIIII with Hooke's 
atom, Sec. II VI with ions, Sec. [V] with a few molecules and 
Sec. I VII with the homogeneous electron liquid. Sec. IVIII 
synthesizes the empirical analysis of Sees. iLTlto lVIl m the 
form of two conjectures. Readers who do not want to go 
through the details of the analysis of different types of 
systems can go right to the conclusions, where all essen- 
tial results are summarized. 

Four appendices deal with issues that are loosely re- 
lated to our main argument, or with by-products of our 
analysis that may be interesting in their own right. Ap- 
pendixfA]motivates and presents a simple analytical fit to 
the atomic data analyzed in Sec. [Ill Appendix [Bl contains 
a comparison of the Lieb-Oxford bound with the Levy- 
Perdew bound, appendix [Cl classifies common (and some 
less common) parametrizations of the electron-liquid cor- 
relation energy with respect to the Lieb-Oxford bound, 
and appendix [D] discusses systems that violate the Lieb- 
Oxford bound. 



II. LIEB-OXFORD BOUND IN ATOMS 



TABLE I: Exchange-correlation energy, LDA exchange en- 
ergy, and their ratio A for light atoms. First three rows (He, 
Be and Ne): Near-exact numerical xc and x energies from 
QMC data, evaluated on exact densitiesji^J^ Rows 4 to 20: 
Precise numerical xc energies obtained^ from CI calculations 
and inversion of the KS equation by the ZMP procedure,— 
and exchange-only LDA energies obtained self-consistently on 
LDA-PW92 densities. For He, Be, and Ne these data are very 
similar to those obtained from QMC, illustrating that the re- 
sulting values of A are robust, and not contaminated with 
inaccuracies due to approximate densities M- All values are in 
atomic (Hartree) units. 





Exc 




A 


He 


1.067 


0.8830 


1.208 


Be 


2.770 


2.321 


1.193 


Ne 


12.48 


11.02 


1.132 


He 


1.068 


0.8617 


1.239 


Li 


1.827 


1.514 


1.207 


Be 


2.772 


2.290 


1.210 


B 


3.870 


3.247 


1.192 


C 


5.210 


4.430 


1.176 


N 


6.780 


5.857 


1.158 





8.430 


7.300 


1.155 


F 


10.320 


8.999 


1.147 


Ne 


12.490 


10.967 


1.139 


Na 


14.440 


12.729 


1.134 


Mg 


16.430 


14.563 


1.128 


Al 


18.530 


16.486 


1.124 


Si 


20.790 


18.544 


1.121 


P 


23.150 


20.743 


1.116 


s 


25.620 


22.950 


1.116 


Cl 


28.190 


25.305 


1.114 


Ar 


31.270 


27.812 


1.124 



In order to calculate A[n] we need the exact LDA 
exchange energy, E X DA , and the exact exchange- 
correlation energy E xc , both on the exact density. For 
a few closed-shell atoms, near-exact E xc values and den- 
sities have been obtained by Umrigar and collaborators 
from QMC calculations^ Of course, near-exact is not 
the same as mathematically exact, but the margin of er- 
ror of these QMC data is much smaller than the effects 
we are after in this work. 

The first three rows of Table U compare the near- 
exact xc energies for He, Be and Ne (from Ref. [H, 
as quoted in Ref. [lrf ). the exact LDA exchange ener- 
gies (obtained by evaluating the LDA functional for ex- 
change on exact QMC densities^) and the resulting ratio 
X[n] = E XC /E^ DA . Two trends immediately leap to the 
eye: (i) the values of A are much smaller than the theo- 
retical upper limit Xlo = 2.275, and (ii) X(Z) decreases 
as a function of atomic number Z. 

To explore these emerging trends for a larger data set, 
the comparison is extended in the other rows of Table|T]to 



other atoms. For these atoms apparently no QMC results 
for E xc and the densities are available. The exchange- 
correlation energies reported in rows 4-20 were extracted 
from Ref. [l7l where they were obtained by numerical in- 
version of the Kohn-Sham (KS) equation on CI densities, 
following the Zhao-Morrison-Parr (ZMP) procedure^ 
The values for E^ DA in rows 4-20 were calculated from 
the exact LDA exchange functional and evaluated at self- 
consistent LDA(PW92) densities^ For He, Be and Ne, 
the resulting values of A can be compared to those ob- 
tained from QMC. Cleary, both sets of data differ slightly, 
but this difference is a small fraction of the difference 
between the observed A values and the theoretical up- 
per bound Xlo- Additionally, we have employed ap- 
proximate xc energies obtained from the B88-LYP GGA 
functional ; 12 i 13 which is highly precise for atoms (and, 
unlike similarly precise nonempirical functionals, such as 
PBE GGA and TPSS meta-GGA, does not make use of 
the Lieb-Oxford bound in its construction). 

Figure [T] illustrates the simple and systematic trend 
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FIG. 1: Value of the energy ratio A, defined in Eq. Ji| as a 
function of atomic number Z for light atoms. Full squares: 
near-exact QMC data. Open circles: CI-ZMP data. Crosses: 
supplementary approximate data obtained from the B88-LYP 
functional. The upper limit Alo is not shown in the figure, 
because on this scale it is too far from the actual values of 
A. The two horizontal lines are estimates of A for Hooke's 
atom and for the Z — > co limit of Helium-like ions, discussed 
in Sec. Mil and Sec. IIVI The continuous line represents the 
analytical expression which is further discussed in Ap- 
pendix [A] 



of A as a function of Z, showing that the Lieb-Oxford 
bound with the originally proposed value Alo = 2.275 
is tightest for small Z, but actually rather generous for 
all atoms, typical values of A being about 50% smaller. 
Note also that the QMC data, which are expected to 
be more precise than the CI-ZMP data, systematically 
predict still smaller values of A. For the Ar atom (Z = 
18), and perhaps also the Be atom (Z = 4), we suspect 
that the ZMP procedure has resulted in a less precise 
correlation energy than for the other atoms, because the 
correlation energy predicted for these atoms by the B88- 
LYP functional is much closer to the extrapolation of the 
QMC and other CI-ZMP data than their CI-ZMP values. 
A simple fit to the B88-LYP data is 



X(Z) = 0.993 



0.313 



(5) 



which is represented by the continuous line in Fig. [T] A 
physical motivation for the form of this fit is given in 
Appendix [A] 

From Table |T] and Fig. Q]we conclude that for atoms A 
follows a simple and systematic trend as a function of Z, 
and always remains far from the upper limit Alo — 2.275, 
approaching approximately half of this value as Z — > 1 
and extrapolating to 1 as Z — > oo. These conclusions 



TABLE II: Exact exchange-correlation and LDA exchange en- 
ergies of Hooke's atom, evaluated on the exact densities, for 
different values of the spring constant k. Data for k = 0.25 
and k — 3.6 x 10 -6 (row one and two) are from Ref. [271 . Al- 
ternative data for k = 0.25 (row three) are from Ref. |23| . and 
the k — * oo limit has been extracted in the way described in 
the main text, from data in Ref. l2ll. 



k 




tpLDA 


A 


3.6 x 10" 6 


-0.0259 


-0.0174 


1.49 


0.25 


-0.5536 


-0.4410 


1.255 


0.25 


-0.555 


-0.4410 


1.26 


k — > oo 






1.17 



are robust with respect to the various different ways of 
obtaining the densities and energies.— Additional com- 
parisons with the Levy-Perdew bound [second inequality 
of Eq. (Q}] are made in Appendix [Bj 



III. LIEB-OXFORD BOUND IN HOOKE'S 
ATOM 

Hooke's atom is a model system in which two electrons 
interact via Coulomb's law but are bound to a harmonic 
potential instead of a 1/r potential. It is frequently used 
in discussing approximate xc functionals and other as- 
pects of DFT, and many exact and numerical results for 
it are known ,2° ^^^^^^ 

The Hamiltonian of Hooke's atom is 



H 



Hookc 



(k) 



» ( V « + v2) + i*(r? + r2) + 1 -5: 
2m 2 i"i — 



r 2 1 

(6) 

It describes two interacting electrons confined in space 
by a harmonic potential of strength k. Since the two 
electrons in this atom interact by the same Coulomb in- 
teraction l/r\2 as in real atoms, this system has the same 
xc functional. The replacement of the nuclear potential 
by a harmonic confinement, on the other hand, greatly 
simplifies the solution of the eigenvalue problem posed 
by Eq. ^ This fact has motivated much work on this 
simple, yet nontrivial, model2&22*&2&&2^± 

For the present purposes we are interested in the ratio 
A[n], which we rewrite as 



AN 



E r 



E c [n] 

E^ DA [n] E^ DA [r 



(7) 



Laufer and Krieger— have shown that for k — > oo the 
LDA exchange energy evaluated on the exact density re- 
covers 85.5% of the exact exchange energy of Hooke's 
atom. Hence, for large k (high curvature), 



E: 



LDA 



E x 



0.855. 



(8) 
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TABLE III: Near-exact exchange-correlation energies and 
LDA exchange energies of ions, evaluated on the exact densi- 
ties, from Ref. |2SI and the resulting value of A. The Z — *• oo 
limit is that of Eq. (|10|) . All energies are in Hartree units. 



1 (H~) -0.422893 -0.337 1.25 

4 (Be 2+ ) -2.320902 -1.957 1.186 

10(Ne 8+ ) -6.073176 -5.173 1.174 

80(Hg 78+ ) -49.824467 -42.699 1.1669 

Z^oo 1.166 



TABLE IV: VMC exchange-correlation energies and LDA ex- 
change energies of a few small molecules=£ and the resulting 
values of A. Last row: DMC data for bulk Si, from Ref. |30| . 



molecule 




TpLDA 


A 


C2H2 


-3.840 


-3.428 


1.120 


C2H4 


-4.606 


-4.110 


1.121 


C2H6 


-5.367 


-4.760 


1.128 


Si 2 


-2.028 


-1.762 


1.151 


Si(bulk) 


-33.23 


-27.66 


1.201 



As k — > 00, the correlation energy of Hooke's atom 
rapidly drops to zero, relative to the exchange energy. 
Hence, we can neglect the second term in Eq. and 
estimate \ Hooke from the inverse of Eq. |(5J), which yields 

X^t e = 1.17, (9) 

The value X^°^ e is shown as a horizontal line in Fig. [1] 

Exact data at some finite values of k, including cor- 
relation, have been presented in Refs. [23| andl27l. These 
data are collected in Table [TTJ The two independent cal- 
culations for k = 0.25 are in excellent agreement, and the 
tendency as a function of k is consistent with our esti- 
mate of the k — ► 00 limit. Clearly, A for Hooke's atom is 
very close to its value for real atoms, and far below the 
limiting value Xlo- 



IV. LIEB-OXFORD BOUND IN IONS FROM 
THE HELIUM ISOELECTRONIC SERIES 

Near-exact numerical data for some representatives of 
the Helium isoelectronic series have been obtained from 
Hylleraas wave functions by Umrigar and Gonze,— and 
are displayed in Table Hill Laufer and Krieger— also con- 
sider the large Z limit of the Helium isoelectronic series, 
for which they find E^ DA /E X = 0.8577. From this we 
obtain, by the same reasoning used for Hooke's atom, 

Aftoc « 1.166. (10) 

The value Xz^*oo is shown as a horizontal line in Fig. [TJ 
and also included in Table HTT1 

The trend of the data in Table HTT1 as a function of Z is 
indeed consistent with A^f+oo ~ 1.166. Interestingly, the 
value of A for the He isoelectronic series is very similar 
to that of Hooke's atom, in particular in the limit of 
very strongly confining external potentials (Z — > 00 and 
k — ► 00, respectively). For all values of Z, including 
negative and positive ions, the resulting values of A are 
much smaller than Xlo- 



V. LIEB-OXFORD BOUND IN MOLECULES 

Exchange-correlation energies for the silicon dimer and 
a few small hydrocarbons have been obtained by Varia- 
tional Monte Carlo (VMC) techniques by Hsing et al^ 
These data, together with the resulting values of A, are 
recorded in Table HVl 

The resulting values of A are quite similar to those 
obtained for atoms. Interestingly, A of the hydrocarbons 
is smaller than that of the C atom, whereas A of the 
silicon dimer is a bit larger than that of the Si atom. Also, 
A for the hydrocarbons slowly grows as a function of the 
number of H atoms. If there is any trend as a function of 
electron number, it is very weak, as is demonstrated by 
the bulk limit, for which included in the Table the DMC 
value for bulk Si, from Ref. l30l . 

Unfortunately, this data set may be too small to draw 
any reliable inferences from such trends. What is beyond 
doubt, however, is that the molecular data predict A val- 
ues that are roughly as far away from the limit Xlo a $ 
previously found for atoms and ions. 



VI. LIEB-OXFORD BOUND IN THE 
ELECTRON LIQUID 

The calculations of the previous sections show that for 
localized atomic, ionic and molecular densities the Lieb- 
Oxford bound is rather generous. The strength of the 
original Licb-Oxford argument, however, rests in the fact 
that is holds for arbitrary densities, and not just for cer- 
tain subsets. In order to extend the investigation to a 
completely different class of densities we thus now turn 
to spatially uniform systems. A priori there is no rea- 
son why one would expect similar values of X[n] to the 
ones found for localized density distributions, although 
the value for bulk Si, mentioned at the end of the previ- 
ous section, strongly suggests so. 

The homogeneous electron liquid is, of course, of 
paramount importance for DFT, as the reference system 
on which the construction of the LDA and many GGAs 
and meta-GGAs are based. It is also of interest in its 
own right as a model for the conduction band of simple 
metals and as a many-body system in which effects of 
the particle-particle interaction can be studied without 



5 



the simultaneous presence of complications due to inho- 
mogeneity in the single-body potential. 31 

The per-particle exchange energy of the homogeneous 
electron liquid is 



= ( — : - 



9tt\ 1/3 1 



(ID 



where r s is the usual electron-liquid parameter related 
to the charge density n via r s = (4t™/3)- 1 / 3 . Below 
we adopt the value Do = 0.4581653, but note that the 
constants in the various parametrizations considered be- 
low are not normally known to this number of significant 
digits. For the electron liquid, and only for the electron 
liquid, the LDA is by construction exact, and Eq. ((4]) 
becomes 



A(r s ) 



e x (r s ) + e c (r s ) 
e x (r s ) 



1 



e c (r s ) 
e x (r s ) ' 



(12) 



The per-particle correlation energy of the electron liq- 
uid is not known in closed form, but the PW92 
parametrization 



8,32 



hi 



e c (r s ) = -2c (l + air s ) x 



2c (P 1 rl /2 + (3 2 r s + /3 3 r 3 s /2 + /? 4 rf) 



(13) 



is the best available fit to the Green's function Monte 
Carlo data of Ref. |M In Eq. ([T3|) en = 0.031091, ft = 
7.5956, ft = 3.5875, ft = 1.6382, ft = 0.49294 and <x x = 
0.21370 are determined such as to reproduce the exactly 
known properties and the QMC data for e c (r s )&£2- 

The unique combination of facts that the LDA becomes 
exact for the electron liquid and that e c is known to very 
high precision, allows us to study A as a continuous func- 
tion of the density parameter r s , instead of at isolated 
densities, as in the previous sections. 

The high-density limit is, rigorously, X(r s — > 0) = 1, 
because \im rg ^ > oe c (r s )/e x (r s ) — 0. To determine the 
low-density limit, recall the leading term of the large-r s 
expansion of the correlation energy^ 



e c {r s 



s do 

00 ) = 

r, 



(14) 



where, according to best estimates^, do — 0.43776. This 
gives the electron-liquid limit for A, 



oo) = 1 



-dp/r s 
-D /r s 



= 1.9555. 



(15) 



This limit is only of formal relevance, as at r s « 65 ± 10 
the electron liquid becomes unstable with respect to the 
Wigner crystal^ and the homogeneous phase ceases to 
be the ground state. Thus, the largest physically possible 
value of the uniform electron liquid is X(r s w 65) = 1.65. 
(This value changes only very little, if the older and pre- 



sumably less accurate estimate 



:s:i 



100 ± 20 for the 
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FIG. 2: Value of the ratio A, denned in Eq. as a function 
of the density parameter r s of the homogeneous electron liq- 
uid, using the PW92 parametrization (|13[l for the correlation 
energy. The four horizontal lines correspond, from top to bot- 
tom, to the Lieb-Oxford^ value Xlo, the Chan-Handy^ value 
Xch, the low-density limit A^^, and the value for bulk Si, 
from Table IIVI Note that actual metals are located at the 
very left margin of the plot, which we have extended to un- 
physically large r s values only to display the large remaining 
distance to Xlo- The filled circle indicates the r s value where 
the transition to the Wigner crystal is expected to occur, 34 
and the uniform electron liquid ceases to be the ground state. 



The quantity do can also be evaluated by employing in 
Eq. (fl2|) available parametrizations of the electron-liquid 
correlation energy. From the PW92 parametrization, as 
specified above, one finds, for example, 



\ PW92 (r s 



ai 4tt / 4 

oo) = 1 H — 

' f3i 3 \9tt 



1/3 



= 1.9462. (16) 



liquid-to-crystal transition is used instead. 



A classification of common (and some less common) 
parametrizations of the electron-liquid correlation energy 
with respect to the Lieb-Oxford bound is presented in 
Appendix [Cl 

We note that the best estimate of do still yields a value 
of A(r s — > oo) that is substantially smaller than the value 
Xlo- Since the value X(r s — > oo) is itself an upper limit 
of A(r s ) at all densities of the electron liquid, this implies 
that in the entire range from r s = to l/r s = the 
Lieb-Oxford bound can be substantially tightened. 

This is illustrated in figure [21 which shows a plot 
of A(r s ) resulting from the PW92 parametrization over 
a wide density interval. Clearly, A(r s ) interpolates 
smoothly between the known limits, and remains far from 
its theoretical upper limit Xlo for any density. Figure [2] 
and the value A(r s — >oo) = 1.9555 show that for any 
possible density of the electron liquid the Lieb-Oxford 



6 



2.25- 

2.00- 

1.75 

1.50- 

1.25- 

1.00- 

0.75- 

0.50- 

0.25- 

0.00 — 



— LO 
-■CH 



: He, Li, 
"Be, B, C, 

N, O, F, 

Na, Ne 

Mg, Ar, Al, Si, 

S, P, CI 




Hooke's atom 



molecules 



solids 



FIG. 3: Graphical summary of all investigated systems. Each panel shows, on the same scale, the Lieb-Oxford bound and its 
refinement by Chan and Handy, together with actual values of A obtained from QMC or CI calculations, or exactly. The first 
panel shows that the entire periodic table fits in a small stretch of A values bounded from above by that of the He atom. Ions, 
small molecules, Hooke's atom, and solids - represented by bulk Si and the electron liquid - occupy almost the same region of 
A space. Values above A = 1.35 are found only for unphysical limits of model Hamiltonians (r s — > oo of the electron liquid, or 
k — > of Hooke's atom) . 



bound falls way above the actual value of A. For metallic 
densities, r s w (2... 6), A is even smaller, not passing 
1.333. 

The A value of bulk Si, A = 1.201, included in Table 
IIV1 falls near the center of the A interval obtained for the 
electron liquid in the metallic density range, although 
bulk Si is not a metal. This again indicates that the 
shape of the density distribution is fairly unimportant 
for the value of A, which never seems to come even near 
the maximum Xlo — 2.275. 



VII. CONCLUSIONS 

The two quantities, X[n] and Xlo, considered in this 
paper have different meanings. X[n], as defined in Eq. Q, 
measures the (inverse) weight of LDA exchange in the full 
exchange-correlation energy of an actual physical system, 
or a model of a physical system. We found here that 
across very different types of systems X[n] does vary, but 
not very strongly. Xlo, as defined in Eq. ([3]), is a math- 
ematical upper limit to X[n], which is the same for all 
nonrelativistic three-dimensional systems with Coulomb 
interactions. 

The analysis of these two quantities, performed in the 
preceding sections, can be summarized as follows: The 
Lieb-Oxford bound provides a lower bound on the xc en- 
ergy of any (nonrelativistic three-dimensional Coulomb) 
system, but nature does not necessarily make use of the 
entire range of permitted values, up to Xlo — 2.275. 
Neutral atoms have A near 1 and approach 1 more closely 
as the atomic number Z is increased. In this limit LDA 



exchange thus captures a larger part of the full xc en- 
ergy. Positive and negative ions, Hooke's atom, small 
molecules and bulk Si all have A values that are close to 
those of isolated atoms, and far from Xlo- The electron 
liquid at metallic densities has A close to 1.3, approach- 
ing 1 and 1.65 in the extreme high and low-density limits, 
respectively. 

Figure [3] contains a graphical summary of our analysis. 
Fine details are not visible on the scale of the figure, but 
the overall impression is very clear, and strongly sug- 
gests that the LO bound is too generous. A tendency 
that systems with a more dilute and spread-out density 
distribution produce larger A values is clearly visible for 
Hooke's atom and the electron liquid. For atoms and 
ions, we find that lighter systems produce larger A val- 
ues. For all investigated systems, A is much smaller than 
Xlo = 2.275. Situations leading to the, relatively, largest 
A values can be classified as follows: 

(i) The largest values of A we found in this investigation 
arise in unphysical low-density limits of model Hamilto- 
nians (the k —y limit of Hooke's atom, or the r s — * oo 
limit of the electron liquid). The largest value we have 
found including such limits is X max = 1-9555, for the 
unattainable low-density limit of the electron liquid. 

(ii) If we exclude physically unrealizable limits of a pa- 
rameter approaching infinity or zero, the largest value we 
have found for physically possible parameters is X max = 
1.65, for the very low-density uniform electron liquid, 
right at its transition to the Wigner crystal. 

(iii) The largest value we have found for any actual 
physical system (atoms, ions, molecules, and solids, but 
excluding Hooke's atom and the electron liquid, which 
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are idealized models) is X max = 1-25, for the H ion. 

This sequence of observations suggests two conjectures. 

Strong conjecture: Tightening the bound for all densi- 
ties. The Lieb-Oxford bound can be substantially tight- 
ened, in the mathematical sense. Independent support 
for this conjecture comes from recent numerical work 
of Chan and Handy^ who indeed found that Xlo can 
be replaced by the slightly smaller value Xch = 2.2149. 
The present investigation suggests, however, that a much 
larger reduction may be achievable. 

It is, of course, conceivable that for certain special den- 
sities X[n] comes arbitrarily close to Xlo ( or Xch), but 
if this is not the case for systems as different as uniform 
electron liquids and isolated atoms and ions, and more- 
over, if the distance from the actual A in these systems 
to Xlo is a large fraction of Xlo itself, it becomes a 
very real possibility that Xlo is not the tightest possible 
limit. Still, if it should turn out to be possible to con- 
struct (perhaps pathological) density distributions that 
require keeping the upper limit at A = Xlo (or Xch), 
this would still be compatible with the following weaker 
conjecture. 

Weak conjecture: Tightening the bound for physical 
densities. For physical densities (arising from realis- 
tic Hamiltonians, excluding unphysical limits) the Lieb- 
Oxford bound can be substantially tightened. In fact, our 
collection of data on different systems suggests a limit of 
Xconjec ~ 1.35, instead of Xlo — 2.275. Interestingly, 
and somewhat suggestively, this empirically found value 
of A makes the prefactor C in E xc [n] > — C J d 3 rn 4 ^ 3 
equal to unity. 

Many existing density functionals, such as the PBE 
GGA 9 and the TPSS meta-GGAiS employ the Lieb- 
Oxford bound, with Xlo = 2.275, as a constraint in 
their construction. It should be most interesting to ex- 
plore if a change of A to a somewhat lower value, either 
universally or for classes of systems, has an impact on 
the performance of these functionals in applications in 
electronic-structure calculations . 
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APPENDIX A: FIT TO THE ATOMIC DATA 

The systematic trend displayed by the atomic data 
in Fig. Q] suggests that their behaviour follows a sim- 
ple law as a function of Z. A reasonable form for this 
law can be guessed as follows: The correct value of 
X(Z — > oo ) is probably exactly one, because asymptoti- 
cally exchange dominates correlation, and LDA exchange 
becomes exact £5. For finite Z, we recall that in gener- 
alized Thomas-Fermi theory 3 -^ all energy contributions 
are expanded powers of Z 1 / 3 , so that the energy ratio 
A should also involve such powers. Since X(Z) decreases 
with increasing Z , the simplest expression consistent with 



TABLE V: Comparison of the Levy-Perdew bound and the 
Lieb-Oxford bound. Values of B are taken from Ref. [5|, where 
they are obtained from a modified PW91 functional, evalu- 
ated at Hartree-Fock densities. The values for A of He and 
Ne are QMC data from Table U For Ar and Kr, Ref. H also 
provides values for B, but for Kr there do not seem to exist 
QMC or CI-ZMP data, and for Ar the CI-ZMP value is so far 
off the curve describing the other atoms (see Fig. [1]) that we 
do not trust it. Hence, we evaluated A for Ar and Kr selfcon- 
sistently from the B88-LYP functional, which is highly precise 
for atoms. To illustrate the error that arises from using an 
approximate functional, we have also included, in row 3, the 
value of A predicted by the B88-LYP functional for Ne. 





A 


B/E^ DA 


He 


1.208 


1.900 


Ne 


1.132 


1.920 


Ne 


1.138 




Ar 


1.111 


1.926 


Kr 


1.079 


1.931 



these expectations is 

\(Z)=A 1 + -^L, (Al) 

where A 2 is a prefactor that cannot be determined by 
such generic considerations, and we allow A\ to deviate 
slightly from unity, because the available data points are 
not exact. 

A simple fit of this expression to the data in Fig.[T]pre- 
dicts Ai = 1.03 and A 2 = 0.231 if only the three QMC 
data points are considered, A\ — 0.986 and A 2 = 0.331 
for the CI-ZMP data (except Ar, which, as explained 
above, may not be properly represented by the CI-ZMP 
data), and A\ = 0.993 and A 2 — 0.313 for the approx- 
imate B88-LYP data, up to Ar. Reasuringly, the fit- 
ted values of A\ are close to the theoretical expectation 
A\ = 1. The third of these fits is shown as continuous 
curve in Fig. [1] 

This simple fit already accounts well for the data in 
Fig. [U suggesting that the proposed Z -1 / 3 behaviour is 
quite realistic. Better fits could, of course, be obtained 
by allowing more terms in the fitting function. However, 
our aim here is not to obtain the best possible fit to the 
data points but to illustrate that for atoms the prefactor 
A has a simple and systematic trend as a function of Z . 

APPENDIX B: LEVY-PERDEW BOUND IN 
ATOMS 

For the He, Ne, Ar and Kr atoms we can compare 
the Lieb-Oxford to the Levy-Perdew bound [second in- 
equality of Eq. (JTJl. Our analysis is based on the data 
in Table 1 of Ref. [a, which reports values of the quan- 
tity B[n] = lim 7 ^o r 7~ 1 E xc [n^] for one of the best avail- 
able GGAs (which at the time of writing of Ref. HI was 
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a slightly modified PW91 3 -), evaluated on tabulated 
Hartree-Fock densities. To the extent that the PW91 
functional and the tabulated Hartree-Fock densities can 
be trusted as approximations to the exact ones, the quan- 
tity B PW91 [n HF ] should, according to Eq. fl}, provide 
a tighter bound on E xc than the Lieb-Oxford inequality, 
as indeed it was found to do£ 

The corresponding values of the ratio 
B PW91 [n HF ]/E^ DA [n HF ] should thus be smaller 
than A^o but still larger than the actual ratio A, defined 
in our Eq. (g]), thus leading to the chain of inequalities 



E , 



JTjLDA 



A < 



D 



JTjLDA 



< A LO - 



(Bl) 



The data in Table [V] show that this expectation is 
bourne out, and that even the Levy-Perdew bound is 
still considerably above the actual (near exact) value of 
A. Interestingly, the tendencies of A and of B/E^ DA as 
functions of Z are opposite, the former decreasing and 
the latter increasing. 



APPENDIX C: 



ELECTRON-GAS CORRELATION 
ENERGY 



TABLE VI: Coefficient do and ratio A for common, and some 
less common, parametrizations of the electron-liquid correla- 



tion energy. 


T^iincinn al 

A Li 1 ± \j 1 vJ 1 L CL± 




\ XT — r*v~i 1 
a ) 


w 

V V 






BR 




1.0631 


f< t 


Pi OO A 71 


1 fill A 

JL.oz.L4 


vBH 


0.56700 


2.2375 


VWN 


0.41433 


1.9043 


PZ81 


0.42681 


1.9316 


PW92 


0.43352 


1.9462 


EHTY 


1.1189 


oo 


elstat.est. 


0.90000 


2.9644 


near-exact 


0.43776 


1.9555 



functional from the first or second class. The BR func- 
tional, which was fitted to data on the He atom^E is a 
special case, as it predicts a value of do that obeys all 
bounds, but is way off the best available value. Hence, 
we do not recommend the use of this expression for ex- 
tended systems. 



We recall from Sec. lVII that the electron liquid displays 
the largest values of A in the low-density limit, r s — > oo, 
where 

A= 1 + ^-1 + ^— = 1.9555. (CI) 

e x _3_ (Stly/ a 

4tt V 4 1 

Many interpolations and parametrizations of e c (r s ) 
have been proposed over the years. In Table I VII we list 
the values of do and X(r s — > oo) predicted by Wigner's 
original interpolation formula (W),— the modification of 
Wigner's expression by Brual and Rothstein (BR);— the 
parametrizations of Gunnarsson and Lundqvist (GL) 40 
and von Barth and Hedin (vBH)^ which are based on 
perturbation theory; those of Vosko, Wilk and Nusair 
(VWN) 4 ^ Perdew and Zunger (PZSl), 4 ^ and Perdew and 
Wang (PW92),- based on Monte Carlo data; a simple 
electrostatic estimate presented in Ref . HH, and the recent 
proposal by Endo et al. (EHTY) ^ which was specifically 
designed for the r s — > limit. 

The entries in Table [VTI fall in three classes. First, the 
electrostatic estimate and the EHTY parametrization vi- 
olate the Lieb-Oxford bound even in its most generous 
universal form, employing Xlo- This is not a surprise, 
considering the crudeness of the electrostatic estimate 
and the fact that the EHTY parametrization was de- 
signed to work well in the r s — > limit, not the r s — > oo 
limit. Second, the von Barth-Hedin parametrization 
and Wigner's interpolation formula obey the Lieb-Oxford 
bound in its universal form, but violate the stricter 
electron-liquid limit, which a proper LDA must also obey. 
Third, all other functionals are consistent also with this 
stricter limit. We do not recommend the use of any xc 



APPENDIX D: SYSTEMS VIOLATING THE 
LIEB-OXFORD BOUND 

Elsewhere in this paper we have repeatedly refered to 
the Lieb-Oxford bound as universal. This use of the con- 
cept of universality is the same commonly employed in 
DFT: a universal quantity (such as the Hohenberg-Kohn 
functional Fjjk) or property (such as the Lieb-Oxford 
bound) is one that is the same for all systems that share 
common kinetic-energy and interaction-energy operators. 
In particular, such quantities or relations are independent 
of the external potentials. Since the Lieb-Oxford bound 
is universal in this sense, we could confront it, above, 
with data on a wide variety of different systems. 

The kinetic-energy operator changes, e.g, in relativis- 
tic quantum mechanics. Even in nonrelativistic quantum 
mechanics it changes if the dimensionality is reduced. We 
stress that the Lieb-Oxford bound was derived for non- 
relativistic three-dimensional systems^ and it is not clear 
if similar results hold in two or one dimensions, or rela- 
tivistically. If similar bounds can be shown to hold, we 
expect that the prefactor Xlo ( or Clo) will be different. 

The interaction-energy operator changes, e.g., when 
DFT is applied to model Hamiltonians i 45 ' 46 i 47 ' 48 An in- 
teresting example is the Hubbard model, where the inter- 
action is local (acting only between electrons at the same 
site) and spin-selective (acting only between electrons of 
opposite spins). 

For fermions, the wave function of the Hubbard model 
is properly antisymmetrized, but this has no conse- 
quences for the energy, i.e., the exchange energy is rig- 
orously zero. The local-density approximation for the 
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Hubbard model, constructed in Ref. l45l and applied, e.g., 
in Refs. l49ll50ll5ll . respects this property. Hence, the Lieb- 
Oxford bound in its Coulomb-interaction form 

E xc [n] > \ LO E% DA [n\ (Dl) 

cannot hold for the Hubbard model, because the right- 



hand side is rigorously zero, whereas the left-hand side is 
known to be nonzero and negative. Similar conclusions 
hold for other model Hamiltonians whose interaction is 
not of Coulomb form. 



* Electronic address: capelle@ifsc.usp.br 

1 W. Kohn, Rev. Mod. Phys. 71, 1253 (1999). 

2 R. M. Dreizler and E. K. U. Gross, Density Functional 
Theory (Springer, Berlin, 1990). 

3 R. G. Parr and W. Yang, Density-Functional Theory of 
Atoms and Molecules (Oxford University Press, Oxford, 
1989). 

4 J. P. Perdew, A. Ruzsinszky, J. Tao, V. N. Staroverov, G. 
E. Scuseria and G. I. Csonka, J. Chem. Phys. 123, 062201 
(2005). 

5 M. Levy and J. P. Perdew, Phys. Rev. A 48, 11638 (1993). 

6 E. H. Lieb and S. Oxford, Int. J. Quantum Chem. 19, 427 
(1981). 

7 A lower bound on C is also known,— but not immediately 
useful for the construction of approximate functionals: to 
obtain a true lower bound on the (negative) functional E xc 
one requires an upper bound for the (positive) quantity 
C. We therefore focus in this paper on this upper bound, 
Clo, which is the value used as input in the construction 
of nonempirical functionals. 

8 J. P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 (1992). 

9 J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett. 
77, 3865 (1996). 

10 J. Tao, J. P. Perdew, V. N. Staroverov and G. E. Scuseria, 
Phys. Rev. Lett. 91, 146401 (2003); J. Chem. Phys. 119, 
12129 (2003); J. Chem. Phys. 120, 6898 (2004); Phys. Rev. 
B 69, 075102 (2004). 

11 J. P. Perdew and Y. Wang, Phys. Rev. B 33, 8800 (1986). 
J. P. Perdew, Phys. Rev. B 33, 8822 (1986). 

12 A. D. Becke, Phys. Rev. A 38, 3098 (1988). 

13 C. Lee, W. Yang and R. G. Parr, Phys. Rev. B 37, 785 
(1988). 

14 G. K.-L. Chan and N. C. Handy, Phys. Rev. A 59, 3075 
(1999). 

15 C. J. Umrigar and X. Gonze, in High Performance Com- 
puting and its Application to the Physical Sciences, Pro- 
ceedings of the Mardi Gras 1993 Conference, eds. D. A. 
Browne et al. (World Scientific, Singapore, 1993). C. Fil- 
ippi, X. Gonze and C. J. Umrigar, in Recent Developments 
and Applications of Modern Density Functional Theory, 
ed. J. Seminario (Elsevier, Amsterdam, 1996). 

16 F. G. Cruz, K.-C. Lam and K. Burke, J. Phys. Chem. A 
102, 4911 (1998). 

17 R. G. Parr and S. K. Ghosh, Phys. Rev. A 51, 3564 (1995). 

18 Q. Zhao, R. C. Morrison and R. G. Parr, Phys. Rev. A 
50, 2138 (1994). 

19 To make sure that our conclusions are stable with respect 
to changes in the employed parametrizations and numeri- 
cal procedures, the ratio A for neutral atoms was also evalu- 
ated in a variety of other ways. These include (i) obtaining 
the exact xc energies not from QMC or ZMP-CI, but by 
subtracting the result of a self-consistent Hartree calcula- 



tion from total energies obtained by full CI,^ (ii) using 
also for He, Be and Ne approximate xc energies and densi- 
ties, obtained from the B88-LY P 12 ' 13 functional, instead of 
the exact ones, (iii) employing the PZ81 and (iv) the VWN 
parametrizations in the self-consistent LDA calculation of 
the exchange energies and densities, instead of PW92, (v) 
employing no parametrization of e c at all, i.e., obtaining 
the exchange energy on exchange-only densities, and (vi) 
using VMC and DMC correlation energies from Ref. |53j in- 
stead of those from CI. (Note that the correlation energies 
of Refs. [k| andl^l are defined with respect to Hartree- Fock 
exchange, not Kohn-Sham exchange, and thus slightly dif- 
ferent from Kohn-Sham correlation energies.) All these 
changes in the procedure lead only to changes in the final 
results that are much smaller than the difference between 
the resulting values of A and the Lieb-Oxford limiting value 
Alo- 

20 N. R. Kestner and O. Sinanoglu, Phys. Rev. 128, 2687 
(1962). 

21 P. M. Laufer and J. B. Krieger, Phys. Rev. A 33, 1480 
(1986). 

22 D. P. O'Neill and P. M. W. Gill, Phys. Rev. A 68, 022505 
(2003). 

23 S. Kais, D. R. Herschbach, N. C. Handy, C. W. Murray 
and G. J. Laming, J. Chem. Phys. 99, 417 (1993). 

24 S. Ivanov, K. Burke and M. Levy, J. Chem. Phys. 110, 
10262 (1999). 

25 M. Taut, Phys. Rev. A 48, 3561 (1993). 

26 W. Zhu and S. B. Trickey, Phys. Rev. A 72, 022501 (2005). 

27 C. Filippi, C. J. Umrigar and M. Taut, 
J. Chem. Phys. 100, 1290 (1994). 

28 C. J. Umrigar and X. Gonze, Phys. Rev. A 50, 3827 
(1994). 

29 C. R. Hsing, M. Y. Chou and T. K. Lee, Phys. Rev. A 74 
032507 (2006). 

30 A. C. Cancio, M. Y. Chou and R. Q. Wood, Phys. Rev. B 
64 115112 (2001). 

31 G. F. Giuliani and G. Vignale, Quantum Theory of the 
Electron Liquid (Cambridge University Press, 2005). 

32 J. P. Perdew and S. Kurth, in A Primer in Density Func- 
tional Theory, eds. C. Fiolhais, F. Nogueira and M. Mar- 
ques (Springer Lecture Notes in Physics Vol. 620, 2003). 
An earlier version appeared in Density Functionals: The- 
ory and Applications, ed. D. Joulbert (Springer Lecture 
Notes in Physics Vol. 500, 1998). 

33 D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 566 
(1980). 

34 G. Ortiz, M. Harris and P. Ballone, Phys. Rev. Lett. 82, 
5317 (1999). 

35 J. P. Perdew, L. A. Constantin, E. Sagvolden and K. 
Burke, Phys. Rev. Lett. 97, 223002 (2006). 

36 L. Spruch, Rev. Mod. Phys. 63, 151 (1991). 



10 



J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson, 
M. R. Pederson, D. J. Singh and C. Fiolhais, Phys. Rev. 
B 46, 6671 (1992). 

E. Wigner, Trans. Faraday Soc. 34, 678 (1938). 

G. Brual, Jr. and S. M. Rothstein, J. Chem. Phys. 69, 1177 

(1978). 

O. Gunnarsson and B. I. Lundqvist, Phys. Rev. B 13, 4274 
(1976). 

U. von Barth and L. Hedin, J. Phys. C 5, 1629 (1972). 
S. H. Vosko, L. Wilk and M. Nusair, Can. J. Phys. 58, 
1200 (1980). 

J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 
T. Endo, M. Horiuchi, Y. Takada and H. Yasuhara, 
Phys. Rev. B 59 7367 (1999). 

N. A. Lima, M. F. Silva, L. N. Oliveira and K. Capelle, 
Phys. Rev. Lett. 90, 146402 (2003). 

K. Capelle, M. Borgh, K. Karkkainen and S.M. Reimann 

submitted (2007), |cond-mat/0702246"l 

V. L. Libero and K. Capelle, Phys. Rev. B 68, 024423 



(2003) . P. E. G. Assis, V. L. Libero and K. Capelle, Phys. 
Rev. B 71, 052402 (2005). 

48 R. J. Magyar and K. Burke, Phys. Rev. A 70, 032508 

(2004) . 

49 N. A. Lima, L. N. Oliveira and K. Capelle, Europhys. Lett. 
60, 601 (2002). 

50 M. F. Silva, N. A. Lima, A. L. Malvezzi and K. Capelle, 
Phys. Rev. B 71, 125130 (2005). 

51 G. Xianlong, M. Polini, M. P. Tosi, V. L. Campo, Jr., K. 
Capelle and M. Rigol, Phys. Rev. B 73, 165120 (2006). 

52 E. R. Davidson, S. A. Hagstrom, S. J. Chakravorty, 
V. M. Umar, and C. Froese Fischer, Phys. Rev. A 44, 7071 
(1991). S. J. Chakravorty, S. R. Gwaltney, E. R. Davidson, 
F. A. Parpia and C. Froese Fischer, Phys. Rev. A 47, 3649 
(1993). 

53 A. Ma, N. D. Drummond, M. D. Towler and R. J. Needs, 
Phys. Rev. E 71, 066704 (2005). 



